Structural instability of two- and three-dimensional pyrochlore spin 

lattices in high magnetic fields 

Oleg Derzhko 1,2,3 and Johannes Richter 2 

institute for Condensed Matter Physics, National Academy of Sciences of Ukraine, 
1 Svientsitskii Street, L'viv-11, 79011, Ukraine 
2 Institut fur Theoretische Physik, Universitat Magdeburg, 
P.O. Box 4120, D-39016 Magdeburg, Germany 
3 National University "Lvivska Politechnika" , 
12 S. Bandera Street, L'viv, 79013, Ukraine 

February 2, 2008 

Abstract 

We consider the quantum spin-s XX Z Heisenberg antiferromagnet on the two- and three-dimensional 
pyrochlore lattices and examine a spin-Peierls mechanism of lowering the total energy by a lattice dis- 
tortion in a high magnetic field. For the exact eigenstates consisting of several independent localized 
magnons in a ferromagnetic environment we show the existence of a spin-Peierls instability by rigorous 
analytical calculations. In addition we report on exact diagonalization data for finite two-dimensional 
pyrochlore lattices up to N — 64 sites. We discuss a scenario of the field-tuned changes in the ground- 
state properties of the frustrated spin lattices due to the coupling between spin and lattice degrees of 
freedom. 

PACS number(s): 75.10.Jm, 75.45. +j 
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1 Introduction. Independent localized-magnon states 

Antifcrromagnetically interacting Heisenberg spins on geometrically frustrated lattices have attracted much 
attention of physicists during last years. Frustrating interactions may lead to a rich phase diagram both 
on the classical and quantum levels |2 |3J 2] . In the presence of an external magnetic field frustrated 
spin systems exhibit a number of specific properties. In particular, plateaus and jumps can be observed in 
the zero-temperature magnetization curve for such models. Another intriguing property is a magnetic-field 
induced spin-Peierls effect in frustrated systems to be discussed below. 

A particular case of geometrically frustrated lattices is the pyrochlore lattice (i.e. a network of corner- 
sharing tetrahedra) which is especially attractive since it has many experimental realizations. The two- 
dimensional (2D) version of the pyrochlore network, also known as the planar pyrochlore, the checkerboard 
lattice, or the square lattice with crossings, has the same local coordination as in the three-dimensional 
(3D) case, however, it is easier for numerical studies. The studies of this simpler model can be helpful in 
understanding of its 3D analog |51 EJ |S] . 

There is a number of experimental and theoretical papers discussing the effects of magnetoelastic cou- 
plings in geometrically frustrated antiferromagnets. Because of frustration, such systems can have a huge 
degeneracy of the ground state which prevents the system from finding a unique ground state as temperature 



decreases. Any perturbation can have a dramatic effect selecting one ground state over another. Different 
perturbations which can lift the ground-state degeneracy have been discussed in the literature, and a magne- 
toelastic coupling is among them. This coupling allows the spin lattice to distort thus relieving the magnetic 
frustration and lowering the total energy of the system. Much attention in recent years has been paid to 
such spin-Peierls-like phase transitions in the pyrochlore antiferromagnets from the experimental as well as 
from the theoretical side. Thus, inelastic magnetic neutron scattering on the 3D pyrochlore antiferromag- 
net ZnCr2C>4 (the spin-3/2 Cr 3+ ions form the 3D pyrochlore lattice) revealed that a lattice distortion can 
lower the energy driving the spin system into an ordered phase [§]. Further probes of a lattice distortion 
in ZnCr2C>4 using neutron scattering and infrared spectroscopy have been reported in Refs. [lUHllj . NMR 
investigation of Y2M02O7 (the spin-1 Mo 4+ ions form the 3D pyrochlore lattice) gives evidence for discrete 
lattice distortions which reduce the energy of the system ^21- The results of a similar /iSR study for this 
compound have been reported recently in Ref. |13j . A lifting of a macroscopic ground-state degeneracy 
of frustrated magnets through a coupling between spin and lattice degrees of freedom in pyrochlore anti- 
ferromagnets was studied theoretically in Refs. ^] 1151 H^j and the relation of these theoretical studies to 
pyrochlore compounds was discussed. More recently, inspired from the experimental data on ZnCi'204 10 
the authors of Ref. J7| have examined a model of lattice-coupled antiferromagnetic classical spins on the 
pyrochlore lattice. They conclude that a particular pattern - a hexagon contraction - may arise owing to 
a lattice distortion which reduces the energy. While the above mentioned studies refer to zero-field case, 
Penc and coworkers |181 HH] have considered the magnetoelastic coupling at finite magnetic fields and have 
shown that the coupling between classical spin and lattice degrees of freedom in a pyrochlore antiferromagnet 
stabilizes a half-magnetization plateau which occurs at intermediate fields. This conclusion agrees with the 
high-field measurements at up to 47 T for the s = 3/2 Heisenberg pyrochlore antiferromagnet CdCr204 |2(J| . 

In what follows, we want to discuss the impact of localized-magnon states |21M22ll25| on the stability of 
the 2D and 3D pyrochlore lattices against distortions. In contrast to the above mentioned zero-field studies 
on the effects of magnetoelastic couplings in the pyrochlore lattice, high magnetic fields are essential in our 
consideration since only in a high magnetic field the independent localized magnons may become relevant 
for the ground-state properties of the considered systems. The present study extends our earlier results 
referring to the square- kagome and kagome lattices |241 125| . 

To be specific, we consider two geometrically frustrated lattices, namely, the 2D pyrochlore lattice (Figs. 
^121) and the 3D pyrochlore lattice (Figs. The Hamiltonian of N quantum spins of length s reads 



Here the sum runs over the bonds (edges) which connect the neighboring sites (vertices) on the spin lattice 
under consideration, J nrn > are the antiferromagnetic exchange constants between the sites n and m, 
A > is the exchange interaction anisotropy parameter, h is the external magnetic field, and S z = ^ n 
is the z-component of the total spin. We assume that all bonds in the lattice without distortion have the 
same length and hence all exchange constants have the same value J. 

We begin with a brief illustration of the independent localized-magnon states in the considered systems 
(the localized magnons in the uniform 2D and 3D pyrochlore lattices were also considered in Refs. [23 and 
|22j. respectively; for a review on localized magnons see Refs. |261 127p . The localized-magnon states are 
excitations above the fully polarized ferromagnetic state, which is the ground state of the antiferromagnet 
in magnetic fields exceeding the saturation field hi- The localized magnons live on restricted areas of the 
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Figure 1: The undistorted two-dimensional pyrochlore lattice. All bonds have the same strength J. In what 
follows we use as shorthands the notations "empty squares" (squares without bond-crossings) and "crossed 
squares" (squares with bond-crossings which are, in fact, projected tetrahedra). 




Figure 2: The distorted two-dimensional pyrochlore lattice. The solid and the dotted bonds have the 
strengths J (l + V2S) and J (l — \/2S) , respectively, the dashed bonds have the strength J (with an accuracy 
up to 5). Note that the lattice volume remains unchanged after such a distortion. 



Figure 3: The {111} slice of the undistorted three-dimensional pyrochlore lattice. All bonds have the same 
strength J. 




Figure 4: (color online). The distorted {111} slice of the three-dimensional pyrochlore lattice. The thick solid 
bonds (thick black) represent the shrinked hexagons and have the strength J (1 + 6), the thin solid and dotted 
bonds (red) have the strength J (l — 5/ (2v / 3)), the dashed bonds (blue) have the strength J (1 — (5/2). The 
normal solid and dotted bonds (normal black) remain as in the undistorted lattice having the strength J. 
Note that the lattice volume remains unchanged after such a distortion. 



lattice, e.g. on the empty squares of the 2D pyrochlore lattice or on the hexagons of 3D pyrochlore lattice. 
They represent exact eigenstates of the frustrated spin lattice and become ground states in high magnetic 
fields. 

In the 2D pyrochlore lattice without distortion the energy of the eigenstate with one magnon localized 
on an empty square is 

E 1 = -2s J + 2s (2s - 1) J A + 4s (4s - 1) JA + (3N - 20) s 2 J A. (2) 

Here the first and the second terms represent the contribution from the bonds along the square, the third 
term represents the contribution from the bonds which connect the square with the environment, and the 
fourth term represents the contribution from the ferromagnetic environment. For n independent localized 
magnons we have instead of 

E n = 3Ns 2 JA-ne 1 , ei=2sJ(l + 3A) (3) 

and n may vary from 1 to n max = N/8 23 . The maximal value of the filling number n max corresponds to 
a dense magnon crystal in which magnons sit on every forth empty square. This magnon crystal state is 
four-fold degenerate and it breaks spontaneously the translational symmetry of the lattice. 

Let us pass to the 3D pyrochlore lattice without distortions. This lattice can be composed from alternating 
kagome-like and triangular planes stacked along the [111] direction. The localized magnons can exist in a 
kagome-like plane (i.e. the {111} slice of the 3D pyrochlore lattice considered in Ref. [2H|) and are trapped on 
hexagons. The {111} slice of the undistorted 3D pyrochlore lattice is shown in Fig. [3J The eigenstate with one 
magnon trapped on a hexagon has the energy given by Eq. @ (however, with the following contributions 
of different bonds: — 2sJ + 2s (3s — 1) JA (along the edges of the hexagon), 4s (6s — 1) JA (along the 
edges connecting the hexagon with the environment), (3iV — 30) s 2 JA (ferromagnetic environment)). For n 
independent localized magnons we again arrive at Eq. @, however, n may vary from 1 to n max = N/12 
|22| . The magnon crystal state corresponds to maximum filling n max . In this magnon crystal state every 
third hexagon in the {111} slice is occupied by magnon producing thus three- fold degeneracy within a slice 
(as for the kagome lattice). Since these three different crystal states can be chosen independently in each 
slice, the degeneracy of the magnon crystal of the whole 3D pyrochlore lattice is much higher and scales at 
least as 3 L where L ~ N 1 / 3 is the number of the {111} slices in the 3D pyrochlore lattice. 

Each magnon decreases S z by 1 and the localized-magnon state with n independent localized magnons 
has S z = Ns — n. These states have the lowest energies in the corresponding sectors of S z |21l I29| . 
Hence, Eq. J3J with n = Ns — S z gives E m i n (S z ) for the 2D and 3D pyrochlore antiferromagnets when 
S 2 = iVs,iVs-l,...,7Vs-7i max . 

In the presence of an external field we must add to the energy J2J| the Zeeman term —hS z = — h (Ns — n). 
Then the ground state belongs to an S z determined by the condition h — dE(S z ) / dS z . At the saturation field 
h\ = t\ the energy of all localized-magnon states with n = 1, . . . , n max is the same and the zero-temperature 
magnetization S z jumps between the saturation value Ns and the value Ns — n max . 

On the basis of general arguments (see also Refs. 0J|2]|S]) we expect that the magnon crystal 

state has gapped excitations and hence the system exhibits a magnetization plateau at S z — Ns — n max 
between hi = hi — Ah and h\. Numerical data support this expectation. Thus, a finite-size analysis of 
the data for the plateau width obtained for the spin- 1/2 isotropic Heisenberg antiferromagnet on finite 
2D pyrochlore lattices with N = 24,32,40,48,64 yields for Ah about 0.2 J. Since we know e.g. from 



the triangular lattice |32l I33| and the square-kagome and kagome lattices [2] EH], that the plateau width 
depends on the spin quantum number s we have also calculated the plateau width Ah for the spin-s isotropic 
Heisenberg system on the 2D pyrochlore lattice with TV = 24 sites for the values of s = 1/2, 1, 3/2, 2, 5/2. 
We find that Aft, scales excellently as ^/s leading to a relative plateau width Ah /hi ~ l/\/s vanishing in the 
classical limit s — > oo. 

Finally, we want to underline that the ground state of both pyrochlore quantum antiferromagnets exhibits 
a huge degeneracy at the saturation field owing to the localized-magnon states. This happens because a 
certain trapping cell of the lattice can be either occupied by a magnon or not whereas the energies of the 
states with different number of magnons are the same. As a result, the degeneracy W, that is the number of 
ways to place magnons on the lattice, grows exponentially with N with a crude estimate W > 2™ max 0] . This 
estimate of the ground-state degeneracy can be improved after mapping the lattice with localized magnons 
onto an auxiliary lattice occupied by hard-core objects [3"H 1551 l3"Hll2"?| . 

2 Lattice instability in high magnetic fields 

We want to examine the lattice stability with respect to a spin-Peierls mechanism in high magnetic fields. 
For that we consider a small lattice deformation which does not violate the conditions for the existence of 
localized-magnon ground states, and analyze the change of the total energy which consists of the magnetic 
and the elastic parts. We search for a maximal gain of the magnetic energy due to lattice deformation. 
For the magnon crystal state the Eqs. (01, @ give a hint how such a deformation may be constructed. 
Actually, for the magnon crystal state the nearest-neighbor spin correlation functions are distributed inho- 
mogeneously. Along the empty square or hexagon hosting a localized magnon the transverse correlations 
( s i 3 i+i) + ( s i s i+i) are negative, but all other nearest-neighbor correlations are positive. Therefore, one may 
expect a maximal gain of magnetic energy by increasing the (antiferromagnetic) bonds on the empty square 
or hexagon which hosts localized magnon and decreasing the (antiferromagnetic) bonds on the attaching 
triangles. The corresponding deformations for the 2D and 3D pyrochlore lattices are discussed below. The 
changes of the magnetic energy due to such deformations are calculated analytically for the values of S z 
corresponding to the magnon crystal. Furthermore we present numerical data for lower values of S z in the 
2D case. 

2.1 Two-dimensional pyrochlore spin lattice 

A pattern of distortions which uses the inhomogeneous distribution of nearest-neighbor spin correlations in 
an optimal way and does not violate the conditions of existence of the magnon crystal state is shown in Fig. 
13 This lattice deformation corresponds to the T2 vibration mode of isolated tetrahedron (see Ref. 
i.e. in every tetrahedron a pair of opposite to each other bonds (belonging to empty squares in Fig. ^) are 
stretched and contracted by the same amount. Introducing the parameter 5 which is proportional to the 
displacement of the sites we find the following changes of the exchange interactions (with an accuracy up to 
linear terms in 8) 



The stronger (weaker) bonds are shown in Fig. [5] by solid (dotted) lines; the bonds having unchanged 
strength are shown in Fig. [21 by dashed lines. As a result, the magnetic energy of the distorted magnon 
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crystal state 

£ IWx (<5) = 3Ns 2 JA - n max (e x + 2y/2s{\ + A)Jtf) (5) 

decreases proportionally to 8. The deformation of the 24 bonds forming one cell belonging to a localized- 
magnon area increases the elastic energy by 87 (s/28) 2 where the parameter 7 is proportional to the elastic 
constant of the lattice. As a result, for the distorted magnon crystal state the elastic energy 167<5 2 n max 
increases proportionally to 8 2 . Combining the magnetic and the elastic energies of the deformed lattices we 
find that the total energy for the magnon crystal state has a minimum at 8 = 8* = (v2/16)s(l + A) J/7 7^ 
that explicitly demonstrates the lattice instability. 

So far our results are rigorous. Now the question arises whether the lattice distortion shown in Fig. 
remains energetically favorable for S z < Ns—n max . To discuss this question we perform exact diagonalization 
for finite lattices using J. Schulenburg's spinpack [H7| . Assuming that the lowest magnetic energy for a given 
value of S z depends on 8 like E m i n (S) = -E m i n (0) + AS P and taking 8 of the order . . . we estimate 
the exponent p numerically. If p is less than 2 (and A < 0) the lattice deformation for the considered finite 
lattice is favorable, otherwise it is not. We have calculated the exponent p for the s — 1/2 isotropic (A = 1) 
Heisenberg system on finite lattices of N — 32, 40 and 48 sites in the two sectors of S* just below the sector 
of the magnon crystal state, i.e. for S z = Ns — n max — 1 and S z = Ns — n max — 2. Note that for N = 48 
with S z = Ns — n max — 2 = 16 the lowest eigenvalues of a Hamiltonian matrix of size 6 • 10 7 x 6 • 10 7 must be 
calculated. We obtain a gain in magnetic energy (i.e. A < 0) in all cases. In the sector S z — hN — n max — 1 
we obtain p = 1.0, 2.0, 2.0 for N = 32, 40, 48, respectively. The "outrider" for N — 32 may be attributed 
to finite-size effects, indeed we find p — 2.0 in the next sectors below S z ~ — rt max — 1. For N = 40 and 
N — 48 the exponent remains p = 2.0 in the sector S z = iiV — n max — 2. 

Thus, our analytical calculations demonstrate a spin-Peierls lattice instability of the 2D pyrochlore anti- 
ferromagnet for the lowest energy states in the sectors of S z = Ns— n max , . . . , Ns— 1. The lattice deformation 
is not favorable in the fully polarized state when S z = Ns. Our numerical calculations give evidence that it 
also disappears when S z < Ns — n max . 

Our findings imply that the 2D pyrochlore quantum antiferromagnet at zero temperature being placed 
in a high magnetic field h > hi becomes distorted when the applied field is lowered to the saturation 
field hi = ei, and the lattice deformation disappears when the field becomes smaller than the field at the 
left endpoint of the plateau /12. It should be noted that a hysteresis phenomenon in the vicinity of the 
saturation field is expected: Starting from a large magnetic field h > hi and then decreasing the field the 
distortion sets in at hi = e x . On the other hand, starting from a field below the saturation field and then 
increasing the field the distorted magnon crystal state survives until the saturation field of the distorted 
lattice hi(8*) = t x + (l/8)s 2 (l + A) 2 J 2 / 7 . 

2.2 Three-dimensional pyrochlore spin lattice 

We proceed along the same lines switching to the 3D pyrochlore lattice. A lattice distortion which does not 
violate the conditions for the existence of localized magnons may be constructed by appropriate deformations 
within the kagome planes of the {111} slices. Such a pattern of distortions is illustrated in Fig. 0] Note that 
the triangular planes and the distance between the triangular and the kagome planes remain unchanged. 
The shrinked hexagons form a x \/3 structure, all other sites of the {111} slice not belonging to shrinked 
hexagons do not change their positions (compare Figs. |3|and0J). The described deformation is similar but 



not identical to the one considered in Ref. The deformation discussed in Ref. ^2] assumes that each 

hexagon contracts towards its center of mass and this distortion does not support the localized magnons. 
Denoting by S the parameter which is proportional to the displacement of the mentioned sites we arrive at 
the following change of the exchange interactions (with an accuracy up to linear terms in S) 



along the edges of the contracted hexagons (thick black lines in Fig. along the out-of- plane triangles 
attached to the contracted hexagons (thin solid and dotted red lines in Fig. 0}, along the in-plane triangles 
attached to the contracted hexagons (dashed blue lines in Fig. 0}, respectively. The exchange interaction 
along the bonds which are not connected to the contracted hexagons (normal solid and dotted black lines in 
Fig. remains unchanged and has the strength J. The magnetic energy of the distorted magnon crystal 
state then reads 



and the corresponding elastic energy is 107<5 2 n max . The total energy achieves a minimum at S — S* — 
(s/20)(2 + (6s + y/Z- 1)A/a/3) J/7 7^ demonstrating the lattice instability. A gain this lattice deformation 
is accompanied by hysteresis phenomenon with the saturation field of the distorted lattice hi(8*) = e± + 



Due to the higher dimensionality the finite-size effects for lattice sizes accessible with exact diagonalization 
are much more important for the 3D pyrochlore lattice and there are no reliable numerical data for the sectors 
below S z = Ns — n max . Therefore, the question about the stability of the lattice against the considered 
deformations for magnetic fields below the plateau belonging to the magnon crystal state remains open for 
the 3D case. 

3 Conclusions 

In the present paper we reported a spin-Pcicrls instability in high magnetic fields of the quantum spin-s 
XXZ Heisenberg antiferromagnet on the 2D and 3D pyrochlore lattices in adiabatic treatment. Since both 
lattices support localized-magnon states which are the lowest energy states in the corresponding sectors 
of S z an instability of the uniform lattice near the saturation field has been proved by rigorous analytical 
calculations. Results for smaller fields have been obtained by exact diagonalization. Although we cannot 
exclude that other distortions (different from the ones shown in Fig. Inland Fig. cannot come into the play, 
we have given arguments that the considered deformations are the most favorable ones taking advantage of 
localized magnons in an optimal way and we have shown explicitly that the undistorted lattice is unfavorable 
in the ground state with S z — Ns — n max , . . . , Ns — 1. 

It is interesting to note that while the lattice deformation for the 2D pyrochlore lattice (Fig. [2J occurs 
simultaneously over all lattice and hence lifts the degeneracy of the ground state at the saturation field, the 
lattice deformation for the 3D lattice (Fig. [3} is composed of local deformations and does not necessarily 
reduce the ground-state degeneracy at the saturation field. 

From experimental point of view, it seems to be important to have a magnetization plateau where the 
deformation occurs, since that gives the possibility to observe the effect in a finite range of h. Observation of 
hysteresis in the magnetization and the deformation of pyrochlore lattices in the vicinity of the saturation field 
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(l/40)s 2 (2 + (6s + V3 - l)A/\/3) 2 J 2 /7- 



can be an experimental sign of the effect discussed theoretically in this paper. For experimental investigations 
of antiferromagnets near saturation field hi one needs materials with comparably small exchange interaction 
J and small value of spin s to be able to reach the high-field regime. Materials with small spin quantum 
number s are more appropriate, also because the plateau near saturation vanishes with increasing s. We 
mention, that in our analysis we dealt with ideal pyrochlore geometry whereas in real compounds deviations 
may be expected. However, there are arguments that the effects of the localized magnons survive under 
slight deviations from the ideal geometry |34| . 

There are now several magnetic materials which can be modeled by the 3D pyrochlore antifcrromagnct 
(for instance, spin-f systems such as ZnV204, MgV204, Y2M02O7, spin-3/2 systems such as ZnCr204, 
CdCi'204, HgCr2C>4; most of them with exchange constants of about 10 2 K) and there is some hope that 
certain features connected with the above presented scenario of a lattice instability driven by a magnetic 
field may be observed in antiferromagnetic compounds with pyrochlore structure. 
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